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Abstract 


The TRMM Multi-satellite Precipitation Analysis (TMPA) provides a calibration-based 
sequential scheme for combining multiple precipitation estimates from satellites, as well as 
gauge analyses where feasible, at fine scales (0.25°x0.25° and 3-hourly). It is available both 
after and in real time, based on calibration by the TRMM Combined Instrument and TRMM 
Microwave Imager precipitation products, respectively. Only the after-real-time product 
incorporates gauge data at the present. The data set covers the latitude band 50°N-S for the 
period 1998 to the delayed present. Early validation results are as follows: The TMPA provides 
reasonable performance at monthly scales, although it is shown to have precipitation rate- 
dependent low bias due to lack of sensitivity to low precipitation rates in one of the input 
products (based on AMSU-B). At finer scales the TMPA is successful at approximately 
reproducing the surface-observation-based histogram of precipitation, as well as reasonably 
detecting large daily events. The TMPA, however, has lower skill in correctly specifying 
moderate and light event amounts on short time intervals, in common with other fine-scale 
estimators. Examples are provided of a flood event and diurnal cycle determination. 
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1. Introduction 


Precipitation displays small-scale variability and highly non-normal statistical behavior that 
requires frequent, closely spaced observations for an adequate representation. Such observations 
are not possible through surface-based measurements over much of the globe, particularly in 
oceanic, remote, and developing regions. Consequently, researchers have come to depend on 
families of sensors flying on a variety of satellites over the last 20+ years for the majority of the 
information used to estimate precipitation on a global basis. While it is possible to create such 
estimates solely from one type of sensor, researchers have increasingly moved to using 
combinations of sensors in an attempt to improve accuracy, coverage, and resolution. The first 
such combinations were performed at relatively coarse scale to ensure reasonable error 
characteristics. For example, the Global Precipitation Climatology Project (GPCP) satellite- 
gauge (SG) combination is computed on a monthly 2.5°x2.5° lat./long. grid (Huffman et al. 
1997; Adler et al. 2003). Subsequently, the scientific community requested that the estimates be 
made available at finer scale, even at the cost of higher uncertainties. Finer scale products 
initiated by the GPCP include the Pentad (Xie et al. 2003) and One-Degree Daily (Huffman et al. 
2001) combination estimates of precipitation. Other research groups have introduced a number 
of fine-scale estimates in the past several years that are now in quasi-operational production (see 
Huffman 2005), including the Climate Prediction Center Morphing algorithm (CMORPH; Joyce 
et al. 2004), the Naval Research Laboratory Global Blended-Statistical Precipitation Analysis 
(NRLgeo; Turk and Miller 2005), the Passive Microwave-calibrated Infrared algorithm (PMIR; 
Kidd et al. 2003), and the Precipitation Estimation from Remotely Sensed Information using 
Artificial Neural Networks (PERSIANN; Sorooshian et al. 2000). 

This paper describes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite 
Precipitation Analysis (TMPA), a new data set that continues the trend toward routine 
computation and distribution of finer-scale precipitation estimates. The primary merged 
microwave-infrared product is computed at the 3-hourly, 0.25°x0.25° lat./long. resolution. In 
common with the GPCP products, the TMPA is designed to combine precipitation estimates 
from various satellite systems, as well as land surface precipitation gauge analyses when 
possible, with the goal that the final product will have a calibration traceable back to the single 
“best” satellite estimate. In the present implementation, the calibration is based on TRMM 
estimates. The TMPA is computed twice as part of the routine processing for TRMM, first as an 
experimental best-effort real-time monitoring product about 6-9 hours after real time, and then as 
a post-real-time research-quality product about 10-15 days after the end of each month. [For 
brevity, these will be referred to as the RT and research products, respectively.] The first has 
been posted to the Web since February 2002, while the second is available from January 1998, 
for a record that totals more than eight years and continues to grow. 

Throughout this paper the research product will be described first, then modifications for 
producing the RT product will be outlined. Section 2 describes the input data sets, while section 
3 describes the algorithms and section 4 gives data set status. Section 5 describes some tests of 
algorithm performance and gives some examples for applying the TMPA to typical applications. 
Section 6 presents concluding remarks. 

2. Input Data Sets 

Most of the coverage in the TMPA depends on input from two different sets of sensors. First, 
precipitation-related passive microwave data are being collected by a variety of low-Earth-orbit 
(LEO) satellites, including the TRMM Microwave Imager (TMI) on TRMM, Special Sensor 
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Microwave/Imager (SSM/I) on Defense Meteorological Satellite Program (DMSP) satellites, 
Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) on 
Aqua, and the Advanced Microwave Sounding Unit B (AMSU-B) on the National Oceanic and 
Atmospheric Administration (NOAA) satellite series. These data have a strong physical 
connection to the hydrometeors that result in surface precipitation, but each individual satellite 
provides a very sparse sampling of the time-space occurrence of precipitation. Even taken 
together, there are significant gaps in the current 3 -hourly coverage by the passive microwave 
estimates. The snapshot in Fig. 1 for a particular 3-hourly period is representative of current 
“full” microwave coverage, averaging about 80% of the Earth’s surface in the latitude band 
50°N-S. The complement of satellite-borne passive microwave sensors has shown steady 
improvement over the span of the TMPA data record (Fig. 2), starting with three satellites in 
1998 that averaged about 40% coverage in each 3-hour period. 

In the current TMPA system, passive microwave pixels from TMI, AMSR-E, and SSM/I are 
converted to precipitation estimates at the TRMM Science Data and Information System 
(TSDIS) with sensor-specific versions of the Goddard Profiling Algorithm (GPROF; Kummerow 
et al. 1996, Olson et al. 1999) for subsequent use in the TMPA. GPROF is a physically based 
algorithm that attempts to reconstruct the observed radiances for each pixel by selecting the 
“best” combination of thousands of numerical-model-generated microwave channel upwelling 
radiances. The associated vertical profiles of hydrometeors then are used to provide an estimated 
surface precipitation rate. The microwave data are screened for contamination by surface 
effects as part of the processing, with marginal contamination denoted as “ambiguous”. 

Passive microwave pixels from AMSU-B are converted to precipitation estimates at the 
National Environmental Satellite Data and Information Service (NESDIS) with operational 
versions of the Zhao and Weng (2002) and Weng et al. (2003) algorithm. Ice Water Path (IWP) 
is computed from the 89- and 150-Ghz channels, with a surface screening that employs ancillary 
data. Precipitation rate is then computed based on the IWP and precipitation rate relations 
derived from cloud model data computed with the NCAR/PSU Mesoscale Model Version 5 
(MM5). The maximum precipitation rate allowed is 30 mm/hr. The AMSU-B algorithm can 
discriminate between precipitating and non-precipitating ice-bearing clouds, but cannot provide 
information on precipitation systems that lack the ice phase. The multi-channel conic-scan 
passive microwave sensors (TMI, AMSR, SSM/I) have a similar limitation over land, so the 
AMSU-B estimates are roughly comparable. But, over ocean the conic scanners also sense 
liquid hydrometeors, providing additional sensitivity, including to precipitation from clouds that 
lack the ice phase. As a result, the AMSU-B estimates over ocean are relatively less capable in 
detecting precipitation over ocean, critically so in subtropical highs. 

The second major data source for the TMPA is the infrared (IR) data that are being collected 
by the international constellation of geosynchronous-Earth-orbit (GEO) satellites. In contrast to 
the sparse temporal sampling of the passive microwave data, the GEO-IR data provide excellent 
time-space coverage. The Climate Prediction Center (CPC) of the National Weather 
Service/NOAA merges the international complement of GEO-IR data into half-hourly 4x4-km- 
equivalent lat./long. grids (hereafter the “CPC merged IR”; Janowiak et al. 2001). The IR 
brightness temperatures (Tb) are corrected for zenith-angle viewing effects and inter- satellite 
calibration differences. 

For research estimates generated prior to the start of the CPC merged IR data set in early 2000 
(see Fig. 2), we use a GPCP data set (also produced at CPC) that contains 24-class histograms of 
GEO-IR Tb data on a 3-hourly, l°xl° lat./long. grid covering the latitude band 40°N-S (hereafter 
the “GPCP IR histograms”; Huffman et al. 2001). This data set also includes grid-box-average 
Geostationary Operational Environmental Satellite (GOES) Precipitation Index (GPI; Arkin and 
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Meisner 1987) estimates computed from LEO-IR data recorded by the NOAA satellite series. 
These LEO-GPI data are used in the TMPA to fill gaps in the GEO-IR coverage, most notably in 
the Indian Ocean sector, where there was no GEO-IR coverage before Meteorological Satellite 5 
(METEOSAT-5) began providing observations there in June 1998. 

All IR-based precipitation estimates share the limitation that the Tb’s primarily sense cloud- 
top temperature, and implicitly cloud height. Arkin and Meisner (1987) showed that such 
information is poorly correlated to precipitation at fine time/space scales, but is relatively well- 
correlated at scales larger than about 1 day and 2.5°x2.5° of lat./long. 

Rounding out this summary of input data sources, the research TMPA also makes use of three 
additional data sources: the TRMM Combined Instrument (TCI) estimate, which employs data 
from both TMI and the TRMM Precipitation Radar (PR), as a source of calibration (TRMM 
product 2B31; Haddad et al. 1997a,b); the GPCP monthly rain gauge analysis developed by the 
Global Precipitation Climatology Centre (GPCC; Rudolf 1993); and the Climate Assessment and 
Monitoring System (CAMS) monthly rain gauge analysis developed by CPC (Xie and Arkin 
1996). 

During the early TMPA design work, the authors realized that they could obtain (restricted) 
access to the requisite microwave and IR data within a few hours of observation time. “Real- 
time” (or more strictly, near-real-time) production makes the estimates useful to several new 
classes of users. At the same time, experience indicated that bias adjustments based on monthly 
gauge data materially improve the accuracy of the estimates (Huffman et al. 1995, among 
others). These considerations led to the two-track approach of computing both RT and research 
products. In the following sections we describe the research product approach, then introduce 
the necessary changes that distinguish the RT algorithm. 

3. TMPA Algorithms 

The TMPA estimates are produced in four stages; (1) the microwave precipitation estimates 
are calibrated and combined, (2) infrared precipitation estimates are created using the calibrated 
microwave precipitation, (3) the microwave and IR estimates are combined, and (4) rain gauge 
data are applied. Figure 3 presents a block diagram of the TMPA estimation procedure. Each 
TMPA precipitation field is best interpreted as the precipitation rate effective at the nominal 
observation time. 

a. Combined microwave estimates 

All of the available passive microwave data are converted to precipitation estimates prior to 
use, then each data set is averaged to the 0.25° spatial grid over the time range ±90 minutes from 
the nominal 3-hourly observation times (00Z, 03Z, . . ., 21Z). All of these estimates are adjusted 
to a “best” estimate using probability matching of precipitation rate histograms assembled from 
coincident data, similar to the probability matched method suggested by Miller (1972) and used, 
for example, by Krajewski and Smith (1991). The research product algorithm takes the TCI as 
the calibrating data source as the result of preliminary validation of Version 6 TRMM products. 
However, the coincidence of TCI with any of the sensors other than TMI is sparse, so we 
establish a TCI-TMI calibration, then apply that to TMI calibrations of the other sensors to 
estimate the TCI-calibrated values. The TCI-TMI relationship is computed on a l°xl° grid for 
each month using that month’s coincident data aggregated on overlapping 3°x3° windows to 
accommodate the somewhat different regional climatologies of the two estimates. 
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Preliminary work showed that the TMI calibrations of the other sensors’ estimates are 
adequately represented by climatologically based coefficients representing large areas. In the 
case of the TMI-SSM/I calibration, separate calibrations are used for five oceanic latitude bands 
(40-30°N, 30-10°N, 10°N-10°S, 10-30°S, 30-40°S) and a single land-area calibration for each of 
the four three-month seasons. The TMI-AMSR-E and TMI-AMSU-B calibrations are set in the 
form of a single climatological adjustment for land and another for ocean. The AMSU-B 
calibration has two additional issues. First, the NESDIS algorithm changed on 31 July 2003, so 
separate sets of calibrations are provided for the two data periods. Second, in both periods the 
AMSU-B fractional occurrence of precipitation in the subtropical highs is notably deficient. 
After extensive preliminary testing, the authors judged it best to develop the ocean calibration in 
regions of significant precipitation and apply it everywhere, recognizing that the resulting fields 
would have a somewhat low bias. In all cases the calibration is a simple match-up of histograms. 
For simplicity, in all cases the calibrations in the 40°-50° latitude belts in both hemispheres are 
taken to be the calibrations that apply just Equator-ward of 40°. 

The calibration interval is chosen to be a month to ensure stability and representativeness, 
except the TMI-AMSR-E calibration is computed using two months for stability. The HQ-IR 
and TCI-TMI calibration interval for the research product is a calendar month, and the resulting 
adjustments are applied to data for the same calendar month. This choice is intended to keep the 
dependent and independent data sets for the calibrations as close as possible in time. 

How do the individual inter-calibrated precipitation estimates compare? All sensor types 
except the AMSR-E recorded views of a storm southeast of Sri Lanka during the 3-hour period 
centered on 00Z of 2 May 2004 (Fig. 4). The enlargement for each subfigure was done without 
interpolations so that individual grid box values are relatively easy to identify. While each 
sensor depicts similar large-scale structure, considerable quantitative disagreement exists at the 
fine scale. Differences in observation times certainly are one factor; the shifts in location of the 
main rain system just south of the Equator are broadly consistent with its westward motion and 
the observation times for each of the sensors. As well, the viewing angles of the LEO satellites 
differ widely for these particular overpasses, giving the possibility that each records a different 
radiometric mix from the same 3 -dimensional scenes. An additional issue for the AMSU-B is 
that its cross-track scanning smooths fields due to the enlargement of pixels at the swath edge. 
Such differences have an important role in creating the large uncertainties at the finest scales, 
discussed below. 

Once the microwave estimates are calibrated for each satellite and audited for >40% 
ambiguous pixels, the grid is populated by the “best” data from all available overpasses. When 
there are multiple overpasses in the 3-hr window for a given grid box, data from TCI, TCI- 
adjusted TMI, TCI-adjusted AMSR-E, and TCI-adjusted SSM/I are averaged together. Tests 
show that the histogram of precipitation rate is somewhat sensitive to the number of overpasses 
averaged together. Accordingly, in the future we expect to test a scheme taking the single “best” 
overpass in the data window period. In contrast, the TCI-adjusted AMSU-B estimates are only 
used if none of the others are available for the grid box due to the detectability deficiency in the 
AMSU-B estimates over ocean (discussed above). Figure 5 (top) provides an example combined 
microwave field. The data voids (gray areas) arise both from the lack of data during the 3 -hour 
period (outside the smoothly arcing swaths of data) and from unfavorable surface conditions 
(northern Japan and southern South America). 

b. Microwave-calibrated IR estimates 
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As noted above, the research product uses two different IR data sets for creating the complete 
record of 3-hrly 0.25° gridded TVs. In the period 1 January 1998 to 7 February 2000, each grid 
box’s histograms in the l°xl° 3-hourly GPCP IR histograms is zenith-angle corrected, averaged 
to a single Tb value for the grid box, and plane-fit interpolated to the 0.25° grid. For the period 
from 7 February 2000 onwards, the CPC Merged IR is averaged to 0.25° resolution and 
combined into hourly files as ±30 minutes from the nominal time. The amount of imagery 
delivered to CPC varies by satellite operator, but international agreements mandate that full 
coverage is provided for the 3-hourly synoptic times (00Z, 03Z, ..., 21Z). Histograms of time- 
space matched HQ precipitation rates and IR T b ‘s, each represented on the same 3-hourly 0.25° 
grid, are accumulated for a month into histograms on a l°xl° grid, aggregated to overlapping 
3°x3° windows, and then used to create spatially varying calibration coefficients that convert IR 
T b ’s to precipitation rates. As in the TCI-TMI calibration for the HQ, the calibration month is 
the calendar month. Screening for ambiguous data described above for the HQ is also applied 
here to the coincident data to control artifacts. 

By design, there is no precipitation when the 0.25°x0.25°-average T b is greater than the local 
threshold value, while increasingly colder T b ’s are assigned increasingly large precipitation rates 
using histogram matching. Those grid boxes that lack coincident data throughout the month, 
usually due to cold-land dropouts or ambiguous editing, are given calibration coefficients by 
smooth-filling histograms of coincident data from surrounding grid boxes. Finally, preliminary 
testing showed that the precipitation rates assigned to the coldest T b ’s by strict probability 
matching tended to show unphysical fluctuations. To ameliorate this effect, a somewhat 
subjectively chosen coldest 0.17% of the T b histogram is specified by a fourth-order polynomial 
fit to a climatology of coldest-0.17%-precipitation rate points around the globe. In each grid box 
a constant is added to the climatological curve to make it is piecewise continuous with the grid 
box’s T b -precipitation rate curve at the 0.17% T b . 

Once computed, the HQ-IR calibration coefficients are applied to each 3-hourly IR data set 
during the month. Figure 5 (upper middle) provides an example microwave-calibrated field; 
spatial coverage by the IR is nearly complete in the latitude band for which we compute 
estimates (50°N-S). Note in Fig. 4 how the IR subfigure qualitatively agrees with the individual 
inter-calibrated microwave estimates, although the area-average precipitation rate and coverage 
by precipitation happen to display higher values in that case. 

c. merged microwave and IR estimates 

The ultimate goal of this project is to provide the “best” estimate of precipitation in each grid 
box at each observation time. It is frequently quite challenging to combine different estimates of 
an intermittent field such as precipitation. The process of combining passive microwave 
estimates is relatively well-behaved because the sensors are quite similar and GPROF is used for 
most retrievals. This is not the case for the microwave and IR fields, as shown in the difference 
field (Figure 5, bottom) for a particular observation time. In addition to the seemingly random 
differences that we would expect, there are also systematic differences (discussed below). 

We apply a simple approach for combining the microwave and IR estimates, namely the 
physically-based combined microwave estimates are taken “as is” where available, and the 
remaining grid boxes are filled with microwave-calibrated IR estimates. This scheme provides 
the “best” local estimate, at the expense of a time series that contains heterogeneous statistics. 
Qualitatively, the example merged field in Figure 5 (lower middle) displays few noticeable data 
boundaries, in part due to the strong spatial variability that real precipitation systems exhibit. 
Note the rain systems associated with then-tropical storm Irene (east of the Bahamas), and 


Huffman et al. TMPA - 3/1/06 


7 


hurricanes Fernando and Greg (southwest and south of Baja California). A mid-latitude low- 
pressure center and trailing front are located north of Hawaii, and late-aftemoon convection is 
decaying in central Africa. 

The difference field (Fig. 5, bottom) facilitates a comparison of the two fields. The 
convective systems such as the African thunderstorms and hurricanes tend to show local 
differences, most likely due to the delay in time between the occurrence of precipitation and the 
growth of cirrus at the top of the storm. The mid-latitude system shows a larger-scale offset, 
which is believed to result from the frontal-scale offset between high-level cirrus ahead of the 
system and precipitation, which is located closer to the low-pressure center and frontal zone. 
Some of the offsets shown in the difference field could also be due to the 3-hour window for 
microwave data. 

Images of the 3-hourly algorithm are best viewed as movie loops, examples of which are 
posted at http://trmm.gsfc.nasa.gov under the button labeled “See More Images, Movies, & 
Accumulation Maps”. 

d. Rescaling to monthly data 

The final step in the research product is the use of rain gauge data. It is highly advantageous 
to include rain gauge data in combination data sets (Huffman et al. 1997, among others). 
However, experience shows that on any time scale shorter than a month the gauge data are not 
reported with sufficient density nor reported with consistent observational intervals to warrant 
direct inclusion in a global algorithm. The authors solved this issue in the GPCP One-Degree 
Daily combination data set by scaling the short-period estimates to sum to a monthly estimate 
that includes monthly gauge data (Huffman et al. 2001). Here, we take a similar approach: All 
available 3-hourly merged estimates are summed over a calendar month to create a monthly 
multi-satellite (MS) product. The MS and gauge are combined as in Huffman et al. (1997) to 
create a post-real-time monthly satellite- gauge combination (SG), which is a TRMM research- 
grade product in its own right (3B43). Then the field of SG/MS ratios is computed on the 0.25° 
grid (with controls) and applied to scale each 3-hourly field in the month, producing the Version 
6 3B42 product. This gauge adjustment step is denoted in Figure 3 with grid-shading. The result 
is to provide the high resolution typical of satellite data and the typically low bias of gauge 
analyses over land. 

e. RT algorithm adjustments 

The RT and research product systems are designed to be as similar as possible to ensure 
consistency between the resulting data sets. One important difference is that the calibrator used 
for the research product, the TCI, is not available in real time. In its absence we use the TMI 
estimates from TRMM. In the future it is planned to test the PR as the standard for calibration 
because it is relatively close to the TCI over land, where the TMI-GPROF estimates are known 
to have several limitations (Kummerow et al. 1996). 

Second, a real-time system cannot reach into the future, so the calibration month is taken as a 
trailing accumulation of 6 pentads + . Before February 2005, the TMI-SSM/I inter-calibration and 
the microwave-IR coefficients are each recomputed at the end of each pentad, while thereafter 


+ A pentad is a 5-day period, except when Leap Day is included in the pentad that encompasses 
it. There are 73 pentads in each year, starting with 1-5 January. 


Huffman et al. TMPA - 3/1/06 


8 


the inter-calibration of individual microwave estimates to the TMI is handled with climatological 
coefficients, and the microwave-IR calibrations are recomputed at each 3-hrly time from the five 
trailing and one current (partial) pentad. The slant-shading on Figure 3 denotes the calibration 
steps in the algorithm block diagram. The change in microwave inter-calibration was made for 
computational simplicity and in preparation for possible loss of TRMM data. The change in 
microwave-IR calibration was intended to better capture rapid changes in the calibration for rare 
heavy rain events that appeared to be poorly handled when the calibration did not include current 
data. 

A third important difference for the RT system is that the monthly gauge adjustment step 
carried out for the research product is not possible. As we gain experience we expect to develop 
a climatological gauge adjustment that can keep the biases of the RT and research products as 
close as possible. 

Figure 6 typifies the relation between the RT and research products. Each data source is 
averaged over the 4-day period 1-4 June 2005, arbitrarily chosen from the epoch with “current” 
algorithm and intended to give a sort of synoptic-scale picture. The spatial averages of the two 
fields are within a few percent of each other, but systematic regional differences are apparent. 
For the most part, these should reflect the differences in the calibrations applied to the two data 
sets. The ITCZ is somewhat lower-amplitude and broader in the post-RT, compared to the RT. 
In the southern extratropics the post-RT is higher almost everywhere, while the picture is 
decidedly mixed in the northern extratropics. Over land the TCI calibration is importantly 
modified by the indirect use of the raingauge analysis in creating the final Version 6 3B42. Low 
values in interior tropical Africa, Central Asia, and the Great Plains in the U.S., as well as higher 
values in Equatorial Amazonia and along the southwestern coast of India are consistent with 
known issues with microwave-based estimates in those regions. 

4. Data Set Status 

The RT system has been running routinely on a best-effort basis in the TRMM Science and 
Data Information System (TSDIS) since late January 2002, and is currently on its third release. 
For simplicity, a fixed latency (currently 9 hours after nominal observation time) has been used 
as the trigger to run the processing. The combined microwave, microwave-calibrated IR, and 
merged microwave-IR estimates, which are labeled 3B-40RT, 3B-41RT, and 3B-42RT, 
respectively, are available from ftp://trmmopen.gsfc.nasa.gov or http://precip.gsfc.nasa.gov. The 
research product system has been developed as the Version 6 algorithm for the TRMM 
operational product 3B-42, although that product only provides the final gauge-adjusted merged 
field. The Version 6 TRMM 3B-43 product has been developed as the post-real-time monthly 
SG described above. Version 6 data are available for January 1998 to the (delayed) present at 
http://lake.nascom.nasa.gov/data/dataset/TRMM/, and an interactive analysis and display 
capability is provided at http://lake.nascom.nasa.gov/Giovanni/tovas. 

5. Comparisons and Examples 

Successive versions of the TMPA-RT have been released to the public since February 2002, 
and various studies have already appeared that examine its performance and utility. For 
example, Katsanos et al. (2004) examined daily-scale performance against gauge data in the 
eastern Mediterranean Sea and found that the TMPA-RT performed best when precipitation was 
relatively frequent, but suffered a bias at high rain rates. In this section we provide a few 
representative comparisons of TMPA performance and show some typical uses. More 
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comprehensive evaluations are in preparation as part of the long-term intercomparison of quasi- 
operational global precipitation data sets over continental- scale regions, the Pilot Evaluation of 
High-Resolution Precipitation Products (PEHRPP), being organized by the International 
Precipitation Working Group (IPWG) of the Coordinating Group for Meteorological Satellites 
(Turk et al. 2006). 

a. Comparison data 

Kwaialein Island. RMI Radar The TRMM Ground Validation (GV) group collects operational 
radar scans from the radar at Kwajalein Island (KWAJ), Republic of the Marshall Islands (RMI), 
for use in TRMM-related studies, to which they apply the recently developed Relative 
Calibration Adjustment (RCA) to correct on-going uncertainty with the radar’s calibration: The 
95 th percentile of each day’s reflectivity histogram collected across known ground clutter 
locations is taken as an index of the radar’s calibration (Marks et al. 2005; Silberstein et al. 
2005). The recalibrated reflectivities are quality-controlled, then matched then against 
coincident raingauge data using the Window Probability Matching Method (Rosenfeld et al. 
1994) to generate Ze-R relationships that vary over time and range from the radar. The resulting 
TRMM 2A-53 GV product is a 2 km x 2 km instantaneous rain map extending 150 km from 
KWAJ (Wolff et al. 2005). For this study the 2-km data were gridded to match the 0.5° 
resolution required for comparison to the buoy gauges and accumulations were generated for 
hourly, 3-hourly, daily, and monthly periods. The 2A-53 data set is nearly unique in its length, 
continuity, level of quality control, and fidelity to the available raingauge data. Even so, users 
need to be aware of residual artifacts, periods of data unavailability, and small-scale fluctuations 
in the Ze-R relationship, all of which introduce uncertainty into the results. 

NWS Melbourne, FL WSR88D Radar Similarly, the TRMM GV group collects operational 
radar scans from the NWS radar at Melbourne, FL (MELB), quality-controls them, then applies 
Window Probability Matching to generate the MELB TRMM 2A-53 GV product. No RCA is 
needed at MELB. 

Australian Gauge Data Daily precipitation accumulations over Australia are provided by the 
Bureau of Meteorology, where reports from up to 6000 gauges are quality-controlled and 
objectively analyzed (Weymouth et al. 1999) to a 0.25° latitude/longitude grid (Weymouth et al. 
1999). The “day” is defined as starting at OOZ. Some portions of interior Australia totally lack 
reporting stations and no analysis is provided in these areas. Elsewhere, the quality of the 
analysis depends on the density of gauge reports (sparser data implies larger errors), the 
complexity of the local terrain (since gauges in mountain areas are generally sited in valleys and 
therefore under-report the true areal average rainfall), and undercatch. The undercatch is 
primarily due to aerodynamic effects at the mouth of the gauge, and so depends on the wind field 
at the time of rainfall. The problem is well-understood (Sevruk 1989), but we lack the requisite 
site and gauge configurations and meteorological data to compute the time-dependent correction 
for each station. The Legates (1987) climatological corrections seem to indicate that the bias 
should be less than 10% in many tropical and warm-season temperate areas. 

Atoll Gauge Data The Comprehensive Pacific Rainfall Database (PACRAIN; Morrissey and 
Green 1991) provides daily and/or monthly island-sited rain gauge data from the western and 
central tropical Pacific Ocean. The “day” is defined locally and carried in the PACRAIN 
metadata. The present study follows the common practice of selecting stations sited on atolls 
and treating them as representative of the open ocean in their region, although this common 
approximation continues to generate controversy. Possible errors include incorrect observation 
times (several of which were identified in this study; Michael Klatt, personal communication), 
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mis-identification of missing values as zeroes, and undercatch. 

Buoy Gauge Data Multi-year observations of rainfall are being recorded on selected buoys in 
the Triangle Trans Ocean Buoy Network (TRITON), Tropical Atmosphere Ocean (TAO), and 
Pilot Research Moored Array in the Tropical Atlantic (PIRATA) programs in the west Pacific, 
east Pacific, and Atlantic Oceans, respectively (see McPhaden et al. 1998). The base observation 
is capacitance-based estimates of water volume in the collection chamber of R. M. Young siphon 
gauges at 1 -minute intervals, which are then converted into filtered 10-minute differences as the 
basic rain rate report (Serra et al. 2001). In the present study the rain rates are averaged over 3- 
hour intervals centered on the synoptic observing times and then all rates below 0.1 mm/h are set 
to zero. This thresholding is necessary to accommodate the noise that broadens the spike in the 
rainrate histogram at zero rain rate into a Gaussian centered slightly below zero; the threshold 
value was chosen in a preliminary examination of the rainrate histogram for several buoy gauges 
with significant rain, but a more careful approach should be taken in the future. As with the 
other gauge data sources, the major source of error is likely a negative bias due to undercatch. 
NERN An integral part of the North American Monsoon Experiment (NAME), which is focused 
on northwestern Mexico, is an array of 87 tipping-bucket gauges that are arranged in transects 
from sea level to the heights of the Sierra Madre known as the NAME Event Rain Gauge 
Network ((NERN; Gochis et al. 2004). The array is specifically laid out to capture the strong 
orographic influence on precipitation that is known to dominate the region’s climatology. 

b. Comparisons 

Monthly Ebert (2005) presented one early set of monthly results that uses the Australian 
raingauge data sets introduced above. In general, she finds that the various satellite estimates, 
including early versions of the TMPA-RT perform best in relatively heavy, convective, warm- 
season regimes, while they perform more poorly in relatively light, cool-season regimes that are 
more mid-latitude in nature. The version of 3B42RT operational in 2004 showed a definite 
seasonal cycle, with Australia-average positive (negative) monthly bias during the warm (cool) 
season and correlation coefficients hovering around 0.50 (0.30) in the warm (cool) season. In 
contrast, we find that the Version 6 3B42 results are significantly closer to the validation data, as 
one might expect with implicit gauge calibration. Calibration coefficients exceed 0.90 for 10 of 
the 12 months in 2004 and there is no obvious seasonal cycle in the bias. As well, 3B42RT 
showed a mean absolute monthly bias of roughly 28% for the entire year 2004, while the Version 
6 TMPA mean absolute monthly bias for 2004 only approaches 9%. The bias over all of 
Australia for 2004 is 6.8%, which is interpreted to mostly reflect the climatological undercatch 
correction that we apply to the TMPA, but which is not present in the Australian analysis. Ebert 
et al. (2006) present further comparisons for Australia, the United States, and the United 
Kingdom that reproduce and extend the results quoted here. 

Turning to the oceans, the atoll and buoy raingauge data sets introduced above may be used as 
the surface reference data because they are not included in processing 3B42/43. The finest 
possible spatial scale for these comparisons is 0.5°x0.5° lat./long. because the buoys are all 
located at the comers of our 0.25° grid boxes. Figure 7 provides the scatter plot for all monthly 
match-ups of (single) raingauge and averaged Version 6 3B43, which is essentially the monthly 
sum of 3B42. The atoll scattergram (Fig. 7a) is quite consistent with the buoy scattergram (Fig. 
7b), showing nearly equitable scatter about the 1:1 line for modest rainrates and a tendency for 
the gauge averages to be higher at higher rainrates. Table 1 lists the corresponding basic 
statistics of average, bias, and RMS. The consistency between the atoll and buoy results tends to 
build confidence that these two rather different data sources are in fact giving a realistic 
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assessment. Further analysis will investigate climate regime dependence, an effect that is 
illustrated by separating the buoy data into locations west and east of the Dateline (Fig. 7b. and 
Table 1). 

In contrast to the point estimates produced by gauges, the estimates from the Kwajalein radar 
are true area averages for a tropical ocean area. The 6-year monthly plot in Fig. 7c shows all the 
0.5° gridbox comparisons within the 100-km range limit for KWAJ, so it only represents one 
area, while the gauge plots cover vast reaches of the tropical Pacific. Nonetheless, the general 
shapes of the scatter plots are similar. The statistics in Table 1 reinforce the conclusion that the 
radar and gauge results are consistent at the monthly scale. Recalling that none of the gauge data 
have been corrected for undercatch, we note that the bias against the atoll gauges is close to that 
against the KWAJ radar estimates, which are calibrated by similar island-based gauges sited 
around Kwajalein Atoll. On the other hand, the bias against the buoy gauges is noticeably 
smaller. We suspect that this is largely due to undercatch. A definitive study of undercatch by 
the buoy gauges is lacking, but Serra et al. (2001) give a range of estimates that is greater than is 
typical for land-based gauges, in part because the buoy gauges are mounted relatively high above 
the surface and in part because the wide-open ocean sites present fewer air-flow obstacles than 
most land sites. This discussion makes it clear that a better evaluation is needed for undercatch 
by the buoy gauges. 

As a comparison using consistent analysis methods, we analyze over-land match-ups to the 
Melbourne WSR-88D (Fig. 7d, Table 1). Although the patterns are similar, the Melbourne land 
results are less biased overall, including less bias for the large number of points at low rain rates. 

Taking an even larger-scale view, consider the monthly time series of precipitation averaged 
over ocean regions for the latitude band 30°N-S (Fig. 8). The Version 6 3B43 (again, essentially 
the sum of the month’s 3B42; shown as a solid line) is closely related to the TCI (gray dash-dot 
line) that is taken as the calibration standard for the passive microwave input fields, except that 
over time the 3B43 shifts to a negative bias of almost 10%. This result seems to stem from the 
AMSU-B rain estimates currently in use. First, as discussed earlier, the AMSU-B algorithm fails 
to detect light rain over oceans, particularly in the subtropical highs. We tempered the influence 
of this AMSU-B deficiency by only using AMSU-B when no other microwave types were 
available, but it is impossible to rescale the AMSU-B rain values when they are all zero. So, the 
successive introduction of AMSU-B satellites during 2000-2002 leads to an increasingly 
negative bias. A second issue for the AMSU-B is that a revised algorithm was applied to data 
observed after 30 July 2003, and the revised algorithm has additional difficulty in identifying 
light precipitation over ocean, even though the overall nature of the rain retrieved is better. The 
authors are presently working to develop a correction approach for the monthly data. A 
comparable issue usually does not arise over land because AMSU-B does not exhibit the same 
problem over land and because the raingauge calibration step keeps the monthly TMPA close to 
the “right” answer. 

The GPCP Version 2 is also included on Fig. 8 for comparison (dashed line), and it tracks 
rather closely with the TCI and Version 6 3B43, except for the drift in 3B43 bias that we ascribe 
to AMSU-B. The calibrator for GPCP over ocean is the Chang-Chiu-Wilheit emission algorithm 
computed for a single-SSM/I time series (Adler et al. 2003), which is quite different than the 
TCI, so the agreement between the two helps build confidence in each result. 

Dailv/Sub-daily Moving to shorter time scales, the TMPA estimates show considerably more 
uncertainty, in common with other short-interval precipitation estimates. It is a challenging area 
of current research to adequately characterize and resolve the uncertainties inherent in such 
estimates. As noted above, the IPWG is organizing an evaluation of high-resolution estimates. 
In this section we focus on illustrating some representative behavior of the TMPA at fine scales. 
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Figure 9 gives a representative time series of the 3-hourly Version 6 3B42 averaged over a 0.5° 
box centered at 5°N,165°E (in the western Pacific ITCZ) for May 2004, as well as the 
TAO/TRITON buoy at that location. Despite the differences in sampling, the two agree on the 
occurrence of most precipitation events. Note that each precipitation event lasts for only a small 
number of 3-hour periods, typical of the convective nature of the ITCZ. The cumulative 
distributions in this particular case are relatively close throughout the month; in other cases a 
different sequence of over- and under-estimates might lead the cumulative distribution to be 
systematically low or high for most of the month. The design goal of the TMPA is that the 
histogram of 3-hourly estimates approximates the histogram of rain that a good validation array 
would observe over the same area and time period. 

Success at this design goal is demonstrated in the matched-data histograms of precipitation at 
the 0.5° 3-hourly scale from Kwajalein for all of 2001 (Fig. 10). The graphs are representative 
of other years of data as well. The important high-end tail matches rather closely, while the 
TMPA estimates show a somewhat more peaked distribution at middle and low values than the 
radar. This is likely the result of the detectability problems noted above in the various satellite 
estimates, particularly from AMSU-B. 

The scattergram for all daily 0.5°x0.5° match-ups between the buoy gauges and V.6 3B42 
estimates for the 7-year period January 1998 to December 2004 is shown in Fig. 1 1. Compared 
to the monthly plots in Fig. 7 there is much more scatter, as is typical of such short-interval 
estimates, particularly when continuous accumulations at a single point (gauge) are scattered 
against sums of area-average snapshots (satellite). For daily rainrates less than about 20 mm d' 1 
the distribution of points is approximately constant along lines perpendicular to the 1:1 line, 
while at higher rainrates the points cluster toward the 1:1 line. It is typical of noisy, calibrated 
estimates that the days with high rates emerge as showing skill first, while the lower values (the 
majority of days) tend to be less skillful as a group. Throughout, the distribution is 
approximately symmetric about the 1 : 1 line, although close inspection reveals that more points 
lie to its right in agreement with the monthly bias results. 

c. Examples 

Among the wide range of possible uses for TMPA estimates, monitoring for extreme 
precipitation events has high societal impact (Negri et al. 2005). In late May 2004 an area of 
disturbed weather moved from south to north across the island of Hispaniola, causing torrential 
rains and flash flooding that caused upwards of 2000 fatalities. This system failed to qualify 
even as a Tropical Depression, but required a significant international relief response. In real 
time the then-operational 3B42RT estimates showed a highly concentrated event in the south- 
central part of the island (Fig. 12, upper left) primarily occurring on 24 May (Fig. 12, upper 
right). Post-analysis seemed to suggest that the estimates were high, and at that time the RT 
scheme only computed a new IR calibration every 5 days, meaning the IR estimates, which 
formed the bulk of the data, were only calibrated by pre-event microwave data. When we shifted 
the IR calibration month by a single 5 -day period to include the flood event, the re-calibrated 
high-end IR estimates were lower by a factor of almost 2 (Fig. 11, lower left). This result 
indicated that extreme events were more sensitive to the calibration than we had expected, and it 
led us to shift the RT system to re-calibrating the IR coeffients every 3 hours. Nonetheless, even 
this new, lower value is symptomatic of a dangerous flood-inducing event. When the Version 6 
3B42 was computed (Fig. 11, lower right), it was similar to the revised RT (Fig. 11, lower left), 
although the research product featured a stronger burst around 10Z on 24 May. Since the 
research product features both a calibration period that encompasses the flood event and much 
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more microwave data than the then-operational 3B42RT, it appears that the IR calibration was 
the dominant source of uncertainty in the RT analysis. 

Another important application of short-interval precipitation estimates such as the TMPA is to 
estimate the diurnal cycle of precipitation. Such information is important for validating model 
performance, as well as providing decision-makers with background data on the most likely 
times for heavy rain events. Comparing the Version 6 3B42 to the NERN in northwestern 
Mexico for the 2004 summer season (June, July, August; Fig. 13), the overall shape of the 
estimated diurnal cycle has good agreement with observations. The main difference is that the 
diurnal maximum through the afternoon and evening is shifted somewhat more than one hour 
later in the TMPA, and is somewhat stronger at its peak. This delay warrants further 
investigation because the plentiful microwave estimates that dominate the research product at 
this time are generally considered closely related to coincident surface precipitation. 

6. Concluding Remarks 

The TMPA has been developed to take advantage of the increasingly rich constellation of 
satellite-borne precipitation-related sensors. Estimates are provided at relatively fine scales 
(0.25°x0.25° 3-hourly) in both post-real and real time to facilitate use by a wide range of 
researchers. However, the errors inherent in the finest scale estimates are large. The most 
successful use of the TMPA data is to take advantage of the fine scales to create averages 
appropriate to the user’s application. 

Looking to future work, the first task for the authors is to further characterize the performance 
of this approach and explore differences between the RT and research products. As part of this 
effort, we plan to participate in PEHRPP (Turk et al. 2006). We also plan to explore 
climatological adjustments to the RT product to minimize its biases against the research version. 
Another important effort is to improve inter-calibration of the microwave-based estimates, 
particularly addressing the oceanic light-precipitation deficiencies in the current AMSU-B input. 

At the same time, we expect to start extending the estimates to the poles by incorporating fully 
global precipitation estimates based on Television Infrared Observation Satellite (TIROS) 
Operational Vertical Sounder (TOVS), Advanced TOVS (ATOVS), and Advanced Infrared 
Sounder (AIRS) data. The best approach to combining the microwave and IR estimates is also a 
topic for future research. It would be helpful to develop a better IR-based algorithm so that the 
combination would not have to reconcile the strong fme-scale differences that currently exist 
between the two. Finally, the study of precipitation in general needs a succinct statistical 
description of how errors in fine-scale precipitation estimates should be aggregated through 
scales up to global/monthly. 

On the instrumentation side there is a concerted effort to provide complete 3-hourly 
microwave data. Most of this effort is focused on the National Aeronautics and Space Admini- 
stration’s Global Precipitation Measurement (GPM) project. Besides simply increasing the 
frequency of coverage, it is planned for GPM to provide a TRMM-like “core” satellite to 
calibrate all of the passive microwave estimates on an on-going basis. We expect the GEO-IR- 
based estimates to have a long-term role in filling the inevitable gaps in microwave coverage, as 
well as in enabling sub-3-hourly precipitation estimates at fine spatial scales. 
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Table Legends 


1. Average, bias, and root-mean-square (RMS) difference statistics comparing monthly 
accumulations reported at single gauges located on atolls and TAO/TRITON buoys, and 
monthly 0.5°x0.5° accumulations of calibrated radar estimates with co-located monthly 
0.5°x0.5° accumulations of the monthly research product (Version 6 3B43) for the 7-year 
period January 1998 to December 2004. All values are expressed as mm/d or percentages, as 
appropriate. 
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Figure Captions 

1. Combined microwave precipitation estimate for the 3-hour period centered on 00Z of 25 
May 2004 in mm/h. Blacked-out areas denote regions that lack reliable estimates, while the 
zero values in the remaining areas are color-coded to depict the coverage by the various 
sensors. The order of precedence for display and corresponding zero color are: TMI (white), 
SSM/I (light grey), AMSR-E (medium grey), and AMSU-B (dark grey). [In the TMPA the 
TMI, SSM/I, and AMSR-E are averaged where overlaps occur.] 

2. Data availability during the TRMM era for satellite sensors used in the TMPA. Solid lines 
denote periods when the data are used in the (Version 6) research product, and dashing 
indicates when they are available, but not used. 

3. Block diagram for both the RT and research product algorithms, showing input data (left 
side), processing (center), output data (right side), data flow (thin arrows) and processing 
control (thick arrows). The items on the slanted shading run asynchronously for the RT 
algorithm, and the items on the grid-shading are only performed for the research product. 
“Best” in the top center shaded box is the TMI GPROF precipitation estimate for the RT 
algorithm and the TMI-PR combined algorithm precipitation estimate for the research 
product. 

4. Single-overpass snapshots from various intercalibrated satellite precipitation estimates for a 
region southeast of Sri Lanka for the 3-hour period centered at 00Z on 2 May 2004. The 
SSM/I and AMSU-B swaths are from the DMSP FI 3 and NOAA-K satellites, respectively. 
For each subimage, the southeastern coast of Sri Lanka is in the upper left comer, the 
Equator is a horizontal dashed line, and longitude 90°E is a vertical dashed line. Black 
denotes missing data, except the latitude/longitude/coast lines are white in data voids and 
black otherwise. The time offset of each image from the nominal 00Z is given in minutes on 
each subimage. 

5. Example of output from the RT algorithm, showing the combined microwave field (HQ, or 
3B40RT; top), microwave-calibrated IR (VAR, or 3B41RT; upper middle), merged 
microwave-IR field (HQVAR, or 3B42RT; lower middle), and the VAR-HQ difference 
(bottom) for the nominal observation time 18Z on 24 September 2002. Grayed-out areas 
lack data for the calculation. 

6. Example accumulations of research (Version 6 3B42; top) and RT (3B42RT; middle) 
products, and the map of differences (bottom) for the period 1-4 June 2005 in mm/d. 

7. Scattergrams comparing monthly research product (Version 6 3B42) precipitation estimates 
to monthly averages reported by single gauges located on a) atolls and b) TAO/TRITON 
buoys for the 7-year period January 1998 to December 2004, as well as to c) Kwajalein radar 
and d) over- land Melbourne radar estimates for the 6-year period January 1999 to December 
2004. All gridding is at 0.5°x0.5° monthly, values are expressed as mm/d, and the buoy 
reports are separated into West and East Pacific (west and east of the Dateline, shown as x’s 
and + ’s, respectively). 

8. Time series of the monthly research product (Version 6 3B43; solid), Version 6 TRMM 
Combined Instrument 3B31 (gray dash-dot), and GPCP Version 2 (dashed) averaged over 
ocean regions in the latitude band 30°N-S. 

9. Time series of the 3-hourly precipitation from the gauge on the TAO/TRITON buoy at 
5°N,165°E (black) and the surrounding 0.5°x0.5° average of the research product (Version 6 
3B42; gray) for May 2004, The 3-hourly data are shown with solid lines (in mm/d) and the 
corresponding cumulative time series are shown with dashed lines (in mm) 
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10. Histograms comparing 3-hourly 0.5 Jx0.5° precipitation rates reported by the Kwajelein radar 
(black) with co-located values from the research product (Version 6 3B42; gray) precipitation 
estimates for 2001. 

1 1 . Scattergrams comparing daily accumulations reported at single gauges located on 
TAO/TRITON buoys with co-located daily 0.5°x0.5° accumulations of the research product 
(Version 6 3B42) precipitation estimates for the 7-year period January 1998 to December 
2004. All values are expressed as mm/d. 

12. Accumulation map (top left) and time series of accumulation averaged within the 250-km 
range ring centered at 19°N,72°W over the period 22-24 May 2004 for the original RT 
(3B42RT; top right), a recalibrated RT (3B42RT; bottom left), and the research product 
(Version 6 3B42; bottom right). The range ring is depicted on the map as the black circle 
roughly centered on Hispaniola. 

13. Comparison of diurnal cycle for the research product (Version 6 3B42; squares) and NAME 
Event Rain Gauge Network (NERRN; dots) for the 2004 summer season (June, July, 
August). A cubic-spline interpolation between research product estimates is shown as a 
dashed line. 
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Table 1. Average, bias, and root-mean-square (RMS) difference statistics comparing monthly 
accumulations reported at single gauges located on atolls and TAO/TRITON buoys, and monthly 
0.5°x0.5° accumulations of calibrated radar estimates with co-located monthly 0.5°x0.5° 
accumulations of the monthly research product (Version 6 3B43) for the 7-year period January 
1998 to December 2004. All values are expressed as mm/d or percentages, as appropriate. 


Gauge Source 

Number of 
comparisons 

Avg. Observed 
Precipitation 

Bias 

(3B43-obs.) 

RMS Difference 

Atolls 

1572 

7.27 

-1.21 (-16.6%) 

3.10 (42.4%) 

All Buoys 

1021 

3.87 

-0.29 (-7.5%) 

1.95 (50.4%) 

W. Pac. Buoys 

323 

7.26 

-0.16 (-2.2%) 

2.34 (32.2%) 

E. Pac. Buoys 

698 

2.32 

-0.35 (-15.2%) 

1.74 (75.0%) 

KWAJ 

504 

5.40 

-0.99 (-18.4%) 

1.94 (36.0%) 

MELB land 

441 

3.26 

+0.23 (+7.1%) 

1.17(35.4%) 
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Figure 1 Combined microwave precipitation estimate for the 3-hour period centered on 00Z of 
25 May 2004 in mm/h. Blacked-out areas denote regions that lack reliable estimates, while the 
zero values in the remaining areas are color-coded to depict the coverage by the various 
sensors. The order of precedence for display and corresponding zero color are: TMI (white), 
SSM/I (light grey), AMSR-E (medium grey), and AMSU-B (dark grey). [In the TMPA the 
TMI, SSM/I, and AMSR-E are averaged where overlaps occur.] 


Huffman et al. TMPA - 2/21/06 


Figs-1 


1 

1998 

2000 

2002 | 

i — ► 

2004 

1 

1 

1 

1 

i 

« 

1 




1 

1 

1 

1 

1 

SSM/I f?13 




SSM/I #14 




1 

SSM/I If 15 




1 

1 

1 

1 

« 

1 

1 

1 

J 

1 

1 

AMSU-B hjl 5 



AMSU-B N16 

1 



AMSU-B Isjl7 

j 


T 

i 

i 

! 

» 

- - - * 
i 

i 

i 



! 

i 

i 

i 


Figure 2 Data availability during the TRMM era for satellite sensors used in the TMPA. Solid 
lines denote periods when the data are used in the (Version 6) research product, and dashing 
indicates when they are available, but not used. 
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Figure 3 Block diagram for both the RT and research product algorithms, showing input data 
(left side), processing (center), output data (right side), data flow (thin arrows) and processing 
control (thick arrows). The items on the slanted shading run asynchronously for the RT 
algorithm, and the items on the grid-shading are only performed for the research product. “Best” 
in the top center shaded box is the TMI GPROF precipitation estimate for the RT algorithm and 
the TMI-PR combined algorithm precipitation estimate for the research product. 
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Figure 4 Single-overpass snapshots from various intercalibrated satellite precipitation estimates 
for a region southeast of Sri Lanka for the 3-hour period centered at OOZ on 2 May 2004. The 
SSM/1 and AMSU-B swaths are from the DMSP F13 and NOAA-K satellites, respectively. For 
each subimage, the southeastern coast of Sri Lanka is in the upper left comer, the Equator is a 
horizontal dashed line, and longitude 90°E is a vertical dashed line. Black denotes missing data, 
except the latitude/longitude/coast lines are white in data voids and black otherwise. The time 
offset of each image from the nominal OOZ is given in minutes on each subimage. 


Huffman et ai. TMPA - 2/21/06 


Figs-4 





Figure 5 Example of output from the RT algorithm, showing the combined microwave field 
(HQ, or 3B40RT; top), microwave-calibrated IR (VAR, or 3B41RT; upper middle), merged 
microwave-IR field (HQVAR, or 3B42RT; lower middle), and the VAR-HQ difference (bottom) 
for the nominal observation time 18Z on 24 September 2002. Grayed-out areas lack data for the 
calculation. 
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Figure 6 Example accumulations of research (Version 6 3B42; top) and RT (3B42RT; middle) 
products, and the map of differences (bottom) for the period 1-4 June 2005 in mm/d. 


Huffman et al. TMPA - 2/21/06 


Figs-6 







Figure 7 Scattergrams comparing monthly research product (Version 6 3B42) precipitation 
estimates to monthly averages reported by single gauges located on a) atolls and b) 
T AO/TRITON buoys for the 7-year period January 1998 to December 2004, as well as to c) 
Kwajalein radar and d) over-land Melbourne radar precipitation estimates for the 6-year period 
January 1999 to December 2004. All gridding is at 0.5°x0.5° monthly, values are expressed as 
mm/d, and the buoy reports are separated into West and East Pacific (west and east of the 
Dateline, shown as x’s and +’s, respectively). 
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Figure 8 Time series of the monthly research product (Version 6 3B43; solid), Version 6 TRMM 
Combined Instrument 3B31 (gray dash-dot), and GPCP Version 2 (dashed) averaged over ocean 
regions in the latitude band 30°N-S. 
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Figure 9 Time series of the 3-hourly precipitation from the gauge on the TAO/TRITON buoy at 
5°N,165°E (black) and the surrounding 0.5°x0.5° average of the research product (Version 6 
3B42; gray) for May 2004, The 3-hourly data are shown with solid lines (in mm/d) and the 
corresponding cumulative time series are shown with dashed lines (in mm). 
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Figure 10 Histograms comparing 3-hourly 0.5Jx0.5° precipitation rates reported by the 
Kwajelein radar (black) with co-located values from the research product (Version 6 3B42; gray) 
precipitation estimates for 2001. 
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Figure 1 1 Scattergrams comparing daily accumulations reported at single gauges located on 
TAO/TRITON buoys with co-located daily 0.5°x0.5° accumulations of the research product 
(Version 6 3B42) precipitation estimates for the 7-year period January 1998 to December 2004. 
All values are expressed as mm/d. 
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Figure 12 Accumulation map (top left) and time series of accumulation averaged within the 250- 
km range ring centered at 19°N,72 0 W over the period 22-24 May 2004 for the original RT 
(3B42RT; top right), a recalibrated RT (3B42RT; bottom left), and the research product (Version 
6 3B42; bottom right). The range ring is depicted on the map as the black circle roughly 
centered on Hispaniola. 
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Figure 13 Comparison of diurnal cycle for the research product (Version 6 3B42; squares) and 
NAME Event Rain Gauge Network (NERRN; dots) for the 2004 summer season (June, July, 
August). A cubic-spline interpolation between research product estimates is shown as a dashed 
line. 
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